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Abstract. 

We present new optical and infrared photometric observa- 
tions and high resolution Ha spectra of the periodic radio star 
LSI+61°303. The optical photometric data set covers the time 
interval 1985-1993 and amounts to about a hundred nights. 
A period of ~26 days is found in the V band. The infrared 
data also present evidence for a similar periodicity, but with 
higher amplitude of variation (0!™2). The spectroscopic obser- 
vations include 16 intermediate and high dispersion spectra 
of LSI+61°303 collected between January 1989 and February 
1993. The Ha emission line profile and its variations are ana- 
lyzed. Several emission line parameters - among them the Ha 
EW and the width of the Ha red hump - change strongly at or 
close to radio maximum, and may exhibit periodic variability. 
We also observe a significant change in the peak separation. 
The Ha profile of LSI+61°303 does not seem peculiar for a 
Be star. However, several of the observed variations of the Ha 
profile can probably be associated with the presence of the 
compact, secondary star. 
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1. Introduction 

The early-type star LSI+61°303 (V 615 Cas) is the optical 
counterpart of the variable radio source GT 0236+610, dis- 
covered during a galactic plane radio survey (Gregory & Tay- 
lor, 1978). Taylor & Gregory (1982) found that this object ex- 
hibits strong radio outbursts with a 26.5 d period. Further ob- 
servations (Taylor & Gregory, 1984) established the currently 
accepted value of 26.496±0.008 d. Typically, radio outbursts 
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peak around phases 0.6-0.8 (Paredes et al., 1990). The spectro- 
scopic radial velocity observations of Hutchings & Crampton 
(1981), hereafter HC81, are in agreement with the radio period, 
and give support to the presence of a companion. In addition, 
they also conclude that the optical spectrum corresponds to a 
rapidly rotating B0 V star, with an equatorial disk and mass 
loss. 

All the radio data available to date on the outburst peak 
flux density provide evidence for a strong modulation, over a 
time scale of 4 yr, in the amplitude of the 26.5 d periodic radio 
outbursts (Gregory et al., 1989; Paredes et al., 1990; Estalella 
et al., 1993). The dependence of radio outbursts flux density on 
frequency, the peak time delay, and the general shape of the ra- 
dio light curves can be modeled as continuous relativistic parti- 
cle injection into an adiabatically expanding synchrotron emit- 
ting source (Paredes et al., 1991). Furthermore, recent VLBI 
observations have provided the first high resolution map of 
LSI+61°303 showing a double sub-arcsec structure (Massi et 
al., 1993). The physical parameters derived from these VLBI 
observations and those of Taylor et al. (1992) are in agreement 
with this model. 

The system was detected as an X-ray source by Bignami et 
al. (1981) and has also been proposed (Gregory & Taylor, 1978; 
Perotti et al., 1980) to be the radio counterpart of the COS B 
7-ray source CG135+01 (Hersem et al., 1977). However, this 
last association is still doubtful due to the large 7-ray error 
box. 

Paredes & Figueras (1986) based on UBVRI photomet- 
ric observations detected optical variability roughly correlated 
with the radio light curve. The amplitude was about 0!™1. A 
model based on deformations of the primary star by a compact 
companion in a eccentric system was initially applied by Pare- 
des (1987) to explain it. Optical variability with time scales of 
days has also been reported by Lipunova (1988) who, in addi- 
tion, found short term nightly fluctuations of some hundredths 
of a magnitude. These short time fluctuations were first ob- 
served by Bartolini et al. (1983). More recently, Mendelson & 
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Mazeh (1989) reported an optical modulation with amplitude 
similar to that found by Paredes & Figueras (1986) and with a 
period of 26.62 ± 0.09 d, near to the radio value, in the John- 
son I band. However, their data set was not sufficient to show 
clearly a similar periodicity at shorter wavelengths. 

The photometric results presented in this paper confirm 
that a ~26 d periodicity is actually present in the V band. In 
addition, the general shape of the visual light curve is very 
similar to that observed by Mendelson & Mazeh (1989) in the 
I band. In the JHK near infrared bands, we find clear evidence 
that a similar periodic modulation is also present, with ampli- 
tude of ~0 m 2. 

On the other hand, intermediate resolution spectroscopic 
data of LSI+61°303 suitable for an analysis of the Ha emis- 
sion line profile are available in literature, at the time of writ- 
ing, only from the early papers of Gregory et al. (1979) and 
HC81. These papers outlined the variations of the Ha profile, 
but were far from reaching any firm conclusion regarding the 
mechanism responsible for such variations. In this work, we 
present 16 spectra (9 of them of high resolution, 0.2-0.44 A 
FWHM) obtained with linear detectors. The new spectra al- 
low a more detailed description of the Ha line profile and of its 
variation. We confirm many of the early findings of Gregory et 
al. (1979) and of HC81 concerning line shifts and Ha EW vari- 
ations. Moreover, we find that other intriguing changes in the 
Ha line profile (noticeably the width of the red hump) occur 
at or close to radio outburst. 

2. Photometric observations and results 

The Johnson photometric observations were made at Calar 
Alto (Almerfa, Spain) with the 1.23 m telescope of the Cen- 
tro Astronomico Hispano-Aleman (CAHA) and the 1.52 m 
telescope of the Observatorio Astronomico Nacional (OAN) 
and at the Observatorio del Roque de los Muchachos (ORM, 
La Palma, Spain), using the 1 m Jacobus Kapteyn telescope 
(JKT). They cover the period 1985-1993 and amount to one 
hundred independent photometric measurements. Both Calar 
Alto telescopes are equipped with a one channel photometer 
with a dry-ice cooled RCA 31034 photomultiplier. The JKT 
observations were made using the People's photometer, with 
two channels, which is equipped with EMI 9658AM photomul- 
tipliers. 

The differential photometry was performed using SAO 
12319 (V=8.79 and 1=7.76), SAO 12327 (V=8.15 and 1=6.90) 
and BD+60°493 (V=8.41 and 1=7.04) as comparison stars. 
The majority of measurements were obtained in the Johnson 
V filter, although some simultaneous I filter observations were 
also taken and will be reported here. Differences of magni- 
tude between comparison stars themselves are constant within 
0!™02. Further details of the observing technique are reported 
in Paredes & Figueras (1986). 

The infrared observations were made at the Teide Ob- 
servatory, (Tenerife, Spain), using the 1.5 m Carlos Sanchez 
telescope (TCS) equipped with the continuously variable filter 
(CVF) photometer. The data were corrected for atmospheric 
extinction and flux-calibrated by comparison with an adequate 
sample of standard stars. 

The results of our Johnson photometric observations are 
given in Table |lj First column indicates the Julian date, second 
and third column are, respectively, the Johnson V and I band 
magnitudes of LSI+61°303. A similar format has been used for 



the infrared observations, whose results are presented in Table 
|i| Figure |l| shows, with the same scale, the optical and infrared 
full data set folded with the radio period of 26.496 d. 

3. Photometric analysis 

3.1. V Photometric periodicity 

In order to try to confirm independently the periodic optical 
modulation reported by Mendelson & Mazeh (1989), a period 
analysis was carried out over this entire data set, amounting to 
105 photometric V measurements over the time interval 1985- 
1993. The period analysis of the data was performed by using 
the phase dispersion minimization (PDM) technique (Stelling- 
werf, 1978). This method consists of assuming a trial period 
and then constructing a phase diagram. The phase interval is 
divided into bins and the variance of the data points is com- 
puted in each bin. The weighted mean of the variances is di- 
vided by the total variance of the data. It can be shown that 
local minima of this function correspond to periods present in 
the data or to multiples of such periods. 

Considering that our minimum sampling rate is about one 
day, our period search was made up to a frequency of 0.5 c d _1 . 
The result of PDM analysis of the V band data is shown in Fig. 
|i| The most significant minimum in the explored frequency 
range (from 0.01 to 0.5 cd~ ) occurs on 0.0387 c d _1 and the 
uncertainty that we associate with this frequency is the fre- 
quency resolution of the complete data set, given by ~ 1/T, 
where T is the total length of the data span and is equal to 
0.0004 c d _1 . This corresponds to a period of 25.8±0.3 d. 

From this analysis, it appears evident that a modulation 
with period of ~26 d is actually present in the Johnson V 
photometry of LSI+61°303. This period is similar to that of 
26.62 d found by Mendelson & Mazeh (1989) in the Johnson 

1 band and to the 26.496 d radio periodicity (Taylor & Gre- 
gory, 1984). In Fig. ^ we have plotted our 105 photometric V 
points folded on the 26.496 d radio period and binned into 10 
bins. Error bars indicate the formal estimate of the uncertainty 
of the mean within each bin. For clarity, the data are plotted 
twice. From this figure we note the similar shape with the op- 
tical light curves presented by Mendelson & Mazeh (1989). In 
particular, the presence of a broad brightness maximun near 
radio active phases 0.5-0.9, and a clear minimum around phase 
0.3 can be appreciated. 

Once the existence of an optical modulation with a pe- 
riod near to 26 d has been established in an independent way 
from Mendelson & Mazeh (1989), it is worth carrying out 
an analysis of all long term photometric data available today 
for LSI+61°303. In this way, we have searched for periodici- 
ties, in the range from 20 to 30 d, the ensemble consisting of 
both our data and the photometric points published by Bar- 
tolini et al. (1983), Lipunova (1988) and Mendelson & Mazeh 
(1989), amounting to 204 nights. Using a frequency step of 

2 x 10~ 7 c d" 1 and a bin structure (5,2), the deepest PDM 
minimum corresponds to 26.5±0.2 d, although it is not very 
prominent. This period value is coincident with the radio pe- 
riod. 

3.2. I Photometry 

Our set of Johnson I band observations, listed in Table 
amount to 43 nights only. With this small data set, it is not 
possible to carry out a feasible periodicity search. In Fig. [I we 
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Table 1. Johnson photometric observations of LSI+61°303 



Julian Date 
(2440000+) 


V 


I 


Julian Day 
(2440000+) 


V 


I 


6264.62 


10.74 


9.26 


7452.68 


10.76 


- 


6265.59 


10.74 


9.25 


7456.44 


10.77 


- 


6267.55 


10.75 


9.24 


7456.63 


10.77 


- 


6268.58 


10.70 


9.20 


7802.60 


10.77 


9.31 


6269.61 


10.72 


9.20 


7803.61 


10.78 


9.32 


6270.62 


10.65 


9.14 


7804.60 


10.79 


9.31 


6271.60 


10.71 


9.20 


7805.59 


10.77 


9.29 


6272.60 


10.70 


9.19 


7806.64 


10.77 


9.29 


6273.61 


10.74 


9.23 


7807.64 


10.77 


9.28 


6274.60 


10.74 


9.24 


7808.62 


10.77 


9.29 


6421.52 


10.77 


9.30 


7809.59 


10.76 


9.27 


6423.34 


10.73 


9.28 


7833.67 


10.75 


9.28 


7056.45 


10.76 


9.27 


7837.59 


10.74 


9.29 


7058.43 


10.74 


9.24 


7838.56 


10.72 


9.28 


7060.58 


10.76 


9.27 


7957.34 


10.74 


9.32 


7063.57 


10.73 


9.24 


8176.43 


10.75 


9.30 


7069.64 


10.72 


9.21 


8180.38 


10.74 


9.28 


7096.66 


10.64 


- 


8238.30 


10.76 


9.31 


7097.56 


10.65 


- 


8274.46 


10.79 


- 


7098.64 


10.70 


- 


8275.49 


10.79 


- 


7125.35 


10.68 


- 


8474.69 


10.79 


9.28 


7128.30 


10.67 


- 


8527.46 


10.75 


9.29 


7147.30 


10.75 


- 


8545.53 


10.78 


- 


7148.30 


10.76 


- 


8620.50 


10.75 


- 


7149.59 


10.71 


- 


8898.65 


10.76 


9.29 


7151.40 


10.73 


- 


9003.46 


10.74 


9.28 


7152.43 


10.72 


- 


9004.47 


10.75 


9.29 


7169.28 


10.74 


- 


9005.45 


10.75 


9.28 


7170.45 


10.70 


- 


7125.35 


10.67 


- 


7169.54 


10.71 


- 


7128.30 


10.68 


- 


7171.56 


10.67 


- 


7168.43 


10.80 


- 


7172.35 


10.71 


- 


7169.30 


10.78 


- 


7172.57 


10.64 


- 


7170.44 


10.78 


- 


7231.34 


10.74 


- 


7367.64 


10.75 


- 


7232.35 


10.70 


- 


7368.64 


10.80 


- 


7381.65 


10.74 


- 


7531.32 


10.73 


- 


7381.65 


10.74 


- 


8092.58 


10.77 


- 


7388.63 


10.79 


- 


8194.49 


10.71 


- 


7389.59 


10.78 


- 


8207.51 


10.76 


- 


7390.59 


10.77 


- 


8208.42 


10.74 


- 


7391.64 


10.77 


- 


8209.35 


10.70 


- 


7392.67 


10.78 


- 


8262.55 


10.75 


9.31 


7393.59 


10.78 


- 


8266.31 


10.73 


9.31 


7394.59 


10.75 


- 


8267.31 


10.74 


9.28 


7395.65 


10.76 


- 


8268.30 


10.74 


9.33 


7420.54 


10.70 




8270.31 


10.76 


9.34 


7421.52 


10.74 




8577.53 


10.79 




7422.56 


10.74 




8610.38 


10.76 




7422.69 


10.74 




8611.41 


10.76 




7423.65 


10.77 




8962.45 


10.79 




7424.64 


10.79 




9034.40 


10.74 




7425.59 


10.79 




9034.40 


10.69 




7426.62 


10.78 











Fig. 1. Optical and infrared observations of LSI+61°303 folded on 
the 26.496 d radio period. Phase zero has been set at JD 2443366.775 
(Taylor & Gregory, 1982). From top to bottom, the V, I, J, H, and K 
band arc plotted. The dots at phase 0.7 are from Elias et al. (1985). 
All observations arc plotted twice. 



4 



Table 2. JHKL' photometric observations of LSI+61°303 



Fig. 2. PDM periodogram of the V band data with bin struc- 
ture (5,5) and frequency step 5 X f0 — 5 c d — 1 . The frequencies are 
in cycles per day. The deepest minimum occurs at a frequency of 
0.0387 c d _1 , corresponding to a period of 25.8 d. 



Fig. 3. Bin average V light curve of LSI+61°303 folded on the 
26.496 d radio period. Phase zero has been set at JD 2443366.775 
(Taylor & Gregory, 1982). Errors bars indicate the formal estimate 
of the uncertainty of the mean within each bin. The continuous line 
is plotted for visual aid. All data is plotted twice. 



show the I band observations folded with the 26.496 d radio 
period. The available data covers the phase interval 0.2-0.9 . 
However, this partial light curve presents the same trends as 
that of V band observations. In particular, a maximum near 
the central radio phases and low emission level near phase 0.2 
is clearly seen. 

3.3. JHK bands photometric periodicity 

Our JHK photometric observations are plotted in Fig. [l] as 
a function of radio phase. Also, we have included two points 
(dots) at phase ~ 0.7 observed by Elias et al. (1985). The 
infrared data available have a good coverage over the full 
radio period, and indicates that the infrared light curves of 
LSI+61°303 also present a modulation similar to that of V 
and I bands. However, the infrared high emission state (orbital 
phases ~0.6-0.9) is broader than in the optical, while the min- 
imum emission state (orbital phases ~0.2-0.4) is deeper and 
narrower. 



Julian Day 
(2440000+) 


J 


H 


K 


L' 


7161.46 


8.76 


8 


39 


8 


06 




7163.45 


8.81 


8 


39 


8 


13 


7.60 


7164.38 


8.77 


8 


38 


8 


15 


7.59 


8279.45 


8.73 


8 


32 


8 


03 




9013.47 


8.62 


8 


25 


7 


88 




8284.51 


8.60 


8 


30 


7 


90 




8491.59 


8.74 


8 


34 


8 


04 




8493.57 


8.80 


8 


35 


8 


05 




8494.60 


8.68 


8 


26 


7 


97 




8496.55 


8.62 


8 


22 


7 


91 




8498.51 


8.57 


8 


18 


7 


87 




8590.35 


8.70 


8 


20 


7 


90 




8592.47 


8.71 


8 


29 


8 


01 




8665.42 


8.57 


8 


16 


7 


89 




8669.46 


8.60 


8 


25 


7 


96 




8856.65 


8.66 


8 


20 


7 


88 




8857.67 


8.62 


8 


21 


7 


92 




8858.65 


8.62 


8 


19 


7 


92 





For a single infrared band, the amount of data accumulated 
by us, 18 nights, could not suffice to carry out a significant pe- 
riod search. In order to overcome this problem, we have merged 
all the JHK photometric points after subtracting their respec- 
tive mean and dividing by the corresponding r.m.s. dispersion 
in each filter. This process is roughly equivalent to having a 
very broad bandpass filter, about 1/xm wide. The two points 
observed by Elias et al. (1985) have also been included. This 
provides a data set of relative normalized infrared magnitudes, 
with 60 measurements over which we have carried a PDM pe- 
riod analysis. The minimum sampling rate is about one day. 
So, the PDM search was carried out up to a frequency of 0.5 
c d _1 . The result of PDM analysis of the merged JHK band 
data is shown in Fig. ^ The most prominent minimum oc- 
curs on 0.0370 ± 0.0003 c d _1 , corresponding to a period of 
27.0 ± 0.3 d. Another nearby deep minimum, with comparable 
significance, is found at 0.0376 ± 0.0003 c d _1 , corresponding 
to a period of 26.6 ± 0.3 d. This implies that an infrared mod- 
ulation, with period similar to the radio period, is also present 
in LSI+61°303. The 60 normalized points used are plotted in 
Fig. |H| as a function of radio phase, computed using the radio 
period value of 26.496 d. 

3.4- Photometric discussion 

The JHK infrared light curves of Fig. |lj showing a deep min- 
imum at phase ~0.3 and a rather flat maximum centered 
around phase ~0.8, are reminiscent of light curves from eclips- 
ing variables. From a rotation velocity value of v sin i ~ 360±25 
km s _1 , HC81 suggest that orbital inclination of LSI+61°303 
is close to 90°. So, this makes the eclipse possibility a rather 
reasonable interpretation. 

In addition, the existence of an IR excess in LSI+61°303 
has been reported by D'Amico et al. (1987) and Elias et al. 
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band, meaning that optical variability should be accounted for 
by a different physical mechanism, probably involving X-ray 
heating of parts of the normal star facing the X-ray source. 



Fig. 4. PDM periodogram of the merged JHK band data with bin 
structure (5,5) and frequency step 5 X 10 — 5 c d — 1 . The frequencies 
are in cycles per day. The deepest minimum occurs at a frequency 
of 0.0370 c d — 1 , corresponding to a period of 27.0 d. 



Fig. 5. JHK infrared observations of LSI+61°303 after substracting 
the mean magnitude of each filter and normalizing with their respec- 
tive r.m.s. dispersion. They are shown folded on the 26.496 d radio 
period. Phase zero has been set at JD 2443366.775. All observations 
are plotted twice. 



(1985). This is, however, a rather common situation in Be 
stars, where the IR excess at micron wavelengths is attributed 
to a dense circumstellar envelope (Slettebak, 1979). This en- 
velope can be also partially responsible for the Ha emission. 
Due to the presence of this envelope, its free-free and free- 
bound opacity is also very likely to absorbe the infrared radi- 
ation from any orbiting companion, thus strongly influencing 
the observed light curve beyond a simple geometrical eclipse. 
A theoretical modelling of the JHK light curves, based on this 
eclipse-attenuation scenario, could yield a determination of the 
system orbital parameters (Marti & Paredes, 1994). 

On the other hand, when considering the optical light 
curves, one sees that the minimum is wider and lasts for about 
half an orbital cycle (see Fig. ^). Then, as suggested by Mendel- 
son & Mazeh (1989), the eclipse explanation at these wave- 
lengths is more difficult to accept. However, due to the ap- 
proximate frequency dependence t v oc f~ 3 , the optical free- 
free and free-bound extinction is actually very low. So, any 
attenuation effects will be difficult to appreciate in the optical 



4. Spectroscopic observations and results 

Sixteen spectra covering the Ha spectral range were collected 
during several observing runs from January 1989 to February 
1993. Table |^ reports the dates of observation and basic infor- 
mation on the instrumental setup. The first column indicates 
the spectrum identification, the second and third the obser- 
vatory and telescope used, the fourth, fifth and sixth contain 
the date, UT time and Julian day of observation, respectively. 
Seventh column is the radio phase. Finally, the eigth and ninth 
columns indicate the dispersion and covered spectral range. All 
spectra were recorded employing CCD detectors. They were 
bias subtracted and flat field corrected using the IRAF or FI- 
GARO package, with the exception of Rozhen spectra which 
have been reduced using the pcIPS software package (Smirnov 
et al., 1992). Only the spectrum obtained on Dec. 27, 1990 
(LP1, see Table |^) at La Palma was flux calibrated. 

In Fig. ^| we show our normalized Ha record of LSI+61°303 
ordered sequentially with radio phase and drawn on the same 
scale with an arbitrary offset. For days with multiple measure- 
ments, only the best profile is shown. The continuum underly- 
ing Ha was rectified to unity employing a spline fitting. This 
normalization procedure is somewhat arbitrary for Rozhen 
Observatory spectra, since the spectral range covered is very 
small, and since most of the Ha wings is probably lost in noise. 

The Ha line profile of LSI+61°303 shows broad wings as 
well as a double peaked core. The Full Width Zero Intensity 
(FWZI) of Ha measured on the Mt. Palomar spectra is « 3100 
km s~ . In addition, in the Asiago spectra, the red hump is 
nearly flat topped and shows a broad shoulder to the red. 

Line parameters were measured on each normalized spec- 
trum and are reported in Table First column gives the spec- 
trum identification. Second, third and fourth columns list the 
heliocentric radial velocity of the blue peak, central dip, and 
red peak, respectively. Fifth and sixth columns are the FWHM 
of the blue and red hump, corrected for instrumental profile, 
while seventh column provides the ratio between the blue and 
red peak intensity above continuum. Eighth column gives the 
total Ha EW and the ninth column lists the EW of the wings. 
Finally, the tenth column contains the EW ratio of the red and 
blue humps. 

The radial velocity, the FWHM and the peak height of the 
B and R humps were measured employing a gaussian fitting. In 
the high resolution spectra, the fitting was done after rebinning 
to a dispersion of ~1 A/pixel. We estimated the contribution 
of the wings by subtracting from the Ha profile a model of the 
core composed of two gaussians. 

5. Spectroscopic analysis 

5.1. Radial velocities 

In Fig. 0, we show the radial velocity difference between the R 
and B peaks as well as the central dip radial velocity, both as a 
function of radio phase. The velocity difference between the Ha 
peaks reaches a maximum close to the radio outburst, during 
an interval of two tenths of radio phase. On the other hand, 
Gregory et al. (1979) noted already that the Balmer central 
dip varied within the range —30/— 80 km s _1 over the period 
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Table 3. Summary of spectroscopic observations of LSI+61°303 



Id. 


Observatory 


Telescope 


Date 




UT 


Julian Day 

(2440000+) 


Radio 
Phase 


Dispersion 

(A/pix) 


Spectral Range 

(A) 


AS1 


Asiago 


1.8m 


1989 Jan 


18 


19"29 m 


7545.3 


0.70 


0.22 


6513-6650 


AS2 


Asiago 


1.8m 


1989 Jan 


19 


lg h 27 m 


7546.3 


0.74 


0.22 


6510-6652 


AS3 


Asiago 


1.8m 


1989 Jan 


19 


20 h 21 m 


7546.3 


0.74 


0.22 


6510-6652 


AS4 


Asiago 


1.8m 


1989 Jan 


21 


19 h 55 m 


7548.3 


0.82 


0.22 


6510-6652 


AS5 


Asiago 


1.8m 


1989 Jan 


21 


20 h 47 m 


7548.4 


0.82 


0.22 


6510-6652 


LP1 


La Palma 


2.5m 


1990 Dec 


27 


01 h U m 


8252.5 


0.40 


0.38 


6492-6672 


LP2 


La Palma 


2.5m 


1991 Jan 


27 


23 h 18 m 


8283.5 


0.57 


0.78 


6400-6800 


LP3 


La Palma 


2.5m 


1991 Aug 


28 


02 h 05 m 


8497.5 


0.64 


0.36 


6472-6680 


MP1 


Mt. Palomar 


1.5m 


1992 Aug 


17 


08 h 05 m 


8851.8 


0.02 


1.00 


6255-6925 


MP2 


Mt. Palomar 


1.5m 


1992 Aug 


17 


os^is'™ 


8851.8 


0.02 


1.00 


6255-6925 


MP3 


Mt. Palomar 


1.5m 


1992 Aug 


18 


09 h 57™ 


8852.9 


0.06 


1.00 


6285-6925 


MP4 


Mt. Palomar 


1.5m 


1992 Aug 


19 


W h 04 m 


8853.9 


0.09 


1.00 


6255-6935 


R01 


Rozhen 


2.0m 


1992 Sep 


03 


23' , 30 m 


8869.5 


0.68 


0.10 


6530-6589 


R02 


Rozhen 


2.0m 


1992 Oct 


09 


2lH5 m 


8905.4 


0.04 


0.10 


6543-6602 


R03 


Rozhen 


2.0m 


1993 Feb 


05 


17 h 4Q m, 


9024.2 


0.52 


0.10 


6533-6594 


R04 


Rozhen 


2.0m 


1993 Feb 


06 


19 h 30 m 


9025.3 


0.56 


0.10 


6534-6594 



Table 4. Observed Ha line parameters of LSI+61°303 



Id. iv (B) iv(dip) iv (R) FWHM(B) FWHM(R) B/R -EW(Ha) EW(wings) EW(R)/EW(B) 
(kms' 1 ) (kms' 1 ) (kms- 1 ) (A) (A) (A) (A) 



AS1 


-227 


-50 


154 


5.10 


6.85 


1.14 


11.46 


1.96 


0.79 


AS2 


-208 


-62 


130 


5.57 


8.67 


1.18 


6.26 


0.61 


0.73 


AS3 


-228 


-56 


125 


5.13 


6.48 


1.25 


6.20 


0.95 


0.89 


AS4 


-219 


-52 


139 


4.64 


7.11 


1.14 


7.14 


0.88 


0.72 


AS5 


-223 


-69 


137 


4.64 


6.90 


1.17 


9.82 


0.96 


0.72 


LP1 


-152 


8 


177 


7.14 


5.64 


0.83 


14.22 


2.21 


0.99 


LP2 


-221 


-64 


110 


5.15 


5.10 


0.74 


15.49 


2.68 


0.71 


LP3 


-230 


-72 


108 


5.20 


5.48 


0.86 


13.09 


3.45 


0.87 


MP1 


- 99 


41 


202 


6.32 


5.52 


0.78 


16.22 


2.32 


0.84 


MP2 


- 96 


35 


205 


5.10 


5.79 


0.83 


17.02 


2.60 


0.74 


MP3 


-161 


-17 


130 


5.49 


5.43 


1.03 


12.69 


1.59 


1.00 


MP4 


-191 


-34 


132 


4.64 


4.91 


1.04 


15.43 


1.43 


0.89 


ROl 


-209 


-69 


110 


4.55 


6.68 


0.79 


18.50 




0.64 


R02 


-193 


-46 


130 


3.65 


5.27 


0.95 


12.13 




0.60 


R03 


-215 


-88 


83 


5.55 


5.52 


0.71 


10.42 




0.54 


R04 


-184 


-76 


105 


5.57 


5.15 


0.77 


11.69 




0.68 



February-March 1978. Our data confirm this variation of radial 
velocity. We also find that the central dip has v r > km s _1 , 
but only in two low resolution spectra. A global shift of the Ha 
line is probably £j 60 km s~\ if we neglect the three "out- 
sider" points that come from low resolution, low significance 
observations. A weighted average of the radial velocities over 
the inverse square of dispersion yields v r (B) « —208 km s _1 , 
iv (dip) w — 60 km s _1 . The variation in the peak separation 
may be due to an increase in v r of the red peak and a decrease 
in the v r of the blue peak. 

5.2. Line width and B/R variability 

In Fig. | we represent the FWHMYR) and FWHM(B) as a 
function of radio phase, while Fig. illustrates the B/R peak 
ratio and EW(Ha) as a function of radio phase also. 

The Ha profile of LSI+61°303 is not peculiar with respect 
to those observed in other Be stars (we can, for instance, com- 



pare LSI+61°303 to the Be stars studied by Slettebak et al. 
(1992)). The EW(Ha) of LSI+61°303 is < 20 A. This value 
sets LSI+61°303 among the Ha-weak Be stars. However, the 
value of the Ha FWZI is, to the best of our knowledge, among 
the largest values observed in Be stars. 

Gregory et al. (1979) described LSI+61°303 as having 
Balmer line emission of variable profile intensity and decre- 
ment. Later, HC81 showed that the radial velocity of the cen- 
tral dip, the emission intensity, and the peak ratio were all 
phase related quantities. The Ha line profile did not change in 
its main features over 14 years. The stability of the Ha pro- 
file is undoubtly remarkable, since Be stars sometimes exhibit 
strong Ha profile changes over timescales of several months. 
Our data confirm that the Ha EW changes strongly in corre- 
spondence of the radio outburst. The minimum value of EW 
(6 A) and the maximum value of peak ratio B/R are observed 
between radio phase 0.7-0.8 (see Fig. 0). 
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Fig. 7. Radial velocity difference, between the R and B peaks, and 
heliocentric radial velocity of the central dip. Both are in km s~*and 
horizontal axis is labeled in radio phase. The size of the dots is 
proportional to the dispersion, i.e., the larger the dot, the higher 
the reliability of the measurement. 



Fig. 6. Normalized Ha line profiles of LSI+61°303 at different ra- 
dio phases given on right-hand side. Horizontal scale is wavelength 
in A and vertical scale is in arbitrary intensity units. The Rozhcn 
Observatory spectra have been convolved with a gaussian of FWHM 
= 0.3 A. Identification codes are revealed in Table 3. 



Fig. 8. FWHM(R) and FWHM(B) as a function of radio phase. 
Black dots represent the FWHM of the fitting gaussians and white 
dots correspond to the hump half maximum width. The size of the 
dots is proportional to the dispersion. 



Fig. 9. B/R, peak ratio and EW(H«) in A as a function of radio 
phase. 



As can be seen in Fig. |8| the FWHM(R) increases from 
a value of ~6 A (~250 km s -1 ), at radio phase ~0.5, to a 
value of ~8 A (~350 km s _1 ) at radio phase ~0.7 . So, the red 
hump seems to become substantially broader near the time 
of radio maximum. On the contrary, the FWHM(B) appears 
to decrease slightly at the same time. In the Asiago spectra 
(obtained close to the radio maximum), the FWHM of the red 
hump is visibly larger by w 100 km s _1 than that of the blue 
one. The exact values of the line width depends somewhat on 
the method employed for the measurement. To ascertain that 
this effect is real we also measured directly the half width of 
the two humps. Although the numbers are somewhat different, 
the same effect is evident in Fig. ^| Also, we have computed 
the flux ratio of the blue and red hump, which happens to be 
,$ 1 at all epochs of observation. This implies that the change 
in peak intensity ratio is mainly due to a variation in width of 
the red hump. 

The flux ratio between line core and line wings is relatively 
constant in all spectra (~ 0.15) but LP3. The Ha profile in LP3 
displays a prominent blue wing, with EW(wings) /EW(core) ~ 
0.26, and asymmetry index of the line wings AI~ —0.37, where 
AI is defined as [EW(Red Wing) -EW(Blue Wing)]/EW(Both 
Wings). In the other low resolution observations, the Ha line 
wings are symmetric within the uncertainties. 

5.3. Spectroscopic discussion 

The most widely accepted explanation of the Ha line pro- 
files in Be stars involves a circumstellar disk-like envelope that 
produces the double-peaked line core. Electron scattered Ha 
photons are expected to produce extended line wings. The 
FWZI of LSI+61°303 Ha is « 3100 km s" 1 > 2wsini ~ 780- 
960 km s _1 (HC81). The excess in the line wings is espe- 
cially evident if we model the core as a sum of two gaussians 
(FWZI(core) ~ 1000 km s" 1 ). If we assume that the line core 
is emitted in a disk, and that the velocity field in the disk is ke- 



plerian, we can estimate the ratio between the inner and outer 
radius of the disk. We obtain R ou t/Rin ~ 9.2, for FWZI(core) 
w 1000 km s _1 . The electron scattering optical depth r es for 
a disk-like geometry can also be computed. The density was 
assumed to depend upon r as n e — n e fi(r / Ri„)~ a , and to fade 
exponentially above and below the symmetry plane of the disk. 
For n e ,o = 10 12 cm -3 , and for a = 2, we find that r es £ 0.3, 
and, if the density n e ,o is ^ 10 12 cm -3 , r es ~ 1. If r e3 is so 
large, and if the temperature of the gas is T e £ 10 4 K, as likely, 
extended line wings of FWZI ~ 3000 km s _1 can be produced 
(Poeckert & Marlborough, 1979). 

The circumstellar disk around the B star should have 
Rout ~ 5 x 10 12 cm, a value similar to the length of the semi- 
major axis of the orbit estimated by HC81. If the eccentricity is 
~ 0.75, the secondary star should cross the circumstellar disk 
and, if the circumstellar disk itself is nearly coplanar to the 
plane of the orbit, even sweep across it while close to perias- 
tron. The resultant accretion could produce the radio outburst. 
A similar model has been proposed for the Be star/X-ray bi- 
nary systems A0538-66 and V0332+52 (e.g., Slettebak, 1988). 
The increase in width of the red hump could be due to a non- 
axisymmetric perturbation in the circumstellar disk, occurring 
close to the outer edge of the disk. For instance, the secondary 
may be crossing the circumstellar shell at that time, close to 
its outer radius (it is interesting to note that the blue hump is 
probably wider than the red one at radio phase « 0.4). How- 
ever, the orbital solution of HC81 suggests that the broaden- 
ing is occurring when the star is close to apoastron, where the 
secondary is probably not in contact with the circumstellar 
disk of the Be star, if the radius of the circumstellar disk is 
Rout ~ 5 x 10 12 cm. 

Alternatively, most of the Ha emission could be associated 
with the compact secondary. An accretion disk around a ~1 
Mq compact object would have R ou t ~ 1 X 10 12 cm, without 
considering broad wings produced by electron scattering. The 
strong decrease in the Ha EW at radio maximum can be ex- 
plained either in terms of obscuration of the disk by the Be star 
(Mendelson & Mazeh, 1989), or in terms of reduced emissivity 
in the disk. 

The JHK light curves point toward an eclipsing binary. 
Hence, the plane of the orbit should be close to the line of 
sight. Since the secondary has probably a mass of only 1/6 
or less the mass of the primary, we expect also a large radial 
velocity oscillation in the peak positions (~ 300 km s _1 ). Our 
radial velocity data are not consistent with this, nor are the 
data obtained by the previous investigators. 

Less clear is the interpretation of the variation of the wings; 
we think that a set of homogeneous observations of sufficiently 
high S/N and resolution are needed to finally reject the pos- 
sibility that the line wings might be emitted by the accretion 
disk of the secondary. 

Even if the accretion disk around the compact star does 
not emit the bulk Ha luminosity, the increase in FWHM(R) 
and in Av r — tv(R) — ?v(B) may be caused by an unresolved 
component whose radial velocity reaches a maximum in cor- 
respondence to the radio maximum. We may expect this if, 
for instance, a cloud of line emitting gas would be ballisti- 
cally ejected along the secondary disk axis. This suggestion is 
appealing, since the radial velocity of this unseen component 
should be £ 200 km s _ at radio maximum, consistent with 
the velocity of the bipolar ejecta of radio plasma partially re- 
solved with VLBI techniques by Massi et al. (1993). A strong 
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analogy could be envisaged with the model proposed by Mar- 
tin & Rees (1979) for SS433. Present data on radial velocity 
and FWHM variation can be explained by a combination of 
global line displacement due to the orbital motion of the Be 
star, and of the v r variation of this unresolved component. The 
red shoulder often present in the Ha profile could be a related, 
higher velocity, feature. 

6. Conclusions 

Our Johnson photometric monitoring of LSI+61°303 has 
shown that this object presents a ~26 d periodic modulation in 
the V band. The shape and amplitude of this modulation are 
similar to those found by Mendelson & Mazeh (1989). In addi- 
tion, the J, H and K observations reported in this paper have 
revealed, for the first time, new evidence of infrared variabil- 
ity with similar trends as seen in the optical, but with higher 
amplitude (0!™2). We have established also that the merged 
JHK light curve exhibits a modulation with period similar to 
the radio period. A possible interpretation of this periodicity 
could involve the eclipse and attenuation of the secondary star 
emission by the Be primary and its envelope. 

It is unclear whether an accretion disk around a compact, 
degenerate companion may contribute to the Ha emission, or 
whether the variations observed are due to perturbations pro- 
duced by the companion on the circumstellar disk of the Be pri- 
mary. This problem persists also because Be stars, as a class, 
are far from being well understood. LSI+61°303 clearly de- 
serves more observations from ground and space. Monitoring 
of the Ha profile at high and intermediate resolution will help 
to solve the main ambiguities left by the present investigation. 
It is also desiderable to obtain a new spectroscopic orbital so- 
lution. 
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